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Abstract 

We classify the complete set of dimension-5 operators relevant for the resonant produc¬ 
tion of a singlet of spin 0 or 2 linearly coupled to the Standard Model (SM). We compute the 
decay width of such states as a function of the effective couplings, and provide the matching 
to various well-motivated New Physics scenarios. We then investigate the possibility that 
one of these neutral resonances be at the origin of the excess in diboson production recently 
reported by the ATLAS collaboration. We perform a shape analysis of the excess under full 
consideration of the systematic uncertainties to extract the width Ttot of the hypothetical 
resonance, finding it to be in the range 26 GeV < Ttot < 144 GeV at 95% C.L. We then 
point out that the three overlapping selections WW, WZ, ZZ reported by ATLAS follow 
a joint trivariate Poisson distribution, which opens the possibility of a thorough likelihood 
analysis of the event rates. The background systematic uncertainties are also included in 
our analysis. We show that the data do not require WZ production and could thus in 
principle be explained by neutral resonances. We then use both the information on the 
width and the cross section, which prove to be highly complementary, to test the effective 
Lagrangians of singlet resonances. Regarding specific models, we find that neither scalars 
coupled via the Higgs-portal nor the Randall-Sundrum (RS) radion can explain the ATLAS 
anomaly. The RS graviton with all matter on the infrared (IR) brane can in principle fit 
the observed excess, while the RS model with matter propagating in the bulk requires the 
presence of IR brane kinetic terms for the gauge fields. 
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1 Introduction 


New particles with TeV masses, neutral under the Standard Model (SM) are a common 
prediction of various New Physics (NP) scenarios. Examples include the Kaluza-Klein 
(KK) graviton and the radion in warped extra dimensions [1], the dilaton in theories of 
strongly coupled electroweak breaking [2], Goldstone bosons of extended composite Higgs 
models [3], mesons and glueballs of strongly-coupled theories [4], extra scalars breaking the 
global symmetry of composite Higgs models [5], Higgs portal models [6], and many more. 
Among the various SM-singlet resonances, those of spin 2 and spin 0 have strikingly similar 
couplings to the SM fields, and it is tempting to treat them in a common framework. 

Recently, the ATLAS collaboration has presented a search for narrow resonances de¬ 
caying to electroweak bosons with hadronic final states using the 8 TeV LHC dataset [7]. 
The weak bosons are highly boosted and are thus reconstructed as a single jet each. A 
moderate but intriguing excess has been observed near the dijet mass rrijj = 2 TeV. It 
is thus an interesting question whether the diboson excess could be explained by neutral 
resonances as those predicted in the above scenarios. 

The goal of this work is thus to present a unified approach for spin-0 and spin-2 
resonances coupled to the SM, and apply it to the search performed in Ref. [7]. In a first 
part, we develop a complete effective field theory (EFT) for neutral resonances of spin 0 
and 2. This general analysis is contained in Sec. 2. As it turns out this EFT consists of only 
few operators, which can further be restricted by theoretically well-motivated assumptions, 
such as approximate flavor and CP conservation. All the different neutral resonances listed 
above then have a simple common description in terms of this effective theory. Explicit 
examples of some of these new physics scenarios are then presented and matched to the 
EFT Lagrangian in Sec. 3. Given the concise description of a large class of models in 
terms of few parameters, our EFT can serve as a model-independent framework that can 
be applied to any search for resonances at the LHC. 

In a second part we then perform a detailed statistical analysis of the ATLAS excess. 
A basic characterisation of the diboson excess is performed in Sec. 4. Local discovery sig¬ 
nificances are computed in both frequentist and Bayesian frameworks, showing a moderate 
evidence for the existence of a signal. The shape of the excess is then analysed, taking 
into account all systematic uncertainties. The total width of the hypothetical resonance 
is found to be 26 GeV < T tot < 144 GeV at 95% C.L. Section 5 contains a comprehen¬ 
sive analysis of the total production rates of the excess. The conditional probabilities for 
tagging a true W, Z and QCD jet as either W or Z are obtained from the ATLAS simula¬ 
tions, and provide the tagging probabilities for the WW, WZ, ZZ selections reported by 
ATLAS. We further observe that these three overlapping selections follow a joint trivariate 
Poisson distribution, which opens the possibility of a thorough likelihood analysis of the 
event rates. The tagging probabilities are checked against the full dataset. The estimation 
of dijet background is treated in a way such that the correlations among the three selections 
are taken into account. The uncertainty on this background estimation is then included as 
a systematic in the total likelihood. Using an actual hypothesis testing, we show that the 
data do not require WZ production and could thus in principle be explained by neutral 
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resonances. 


Finally, in a third part, Sec. 6, we test the effective Lagrangians of neutral resonance 
using both the information from the width and from the cross sections. It turns out that 
these pieces of information imply stringent contraints on the EFT parameter space once 
put together, even after including the uncertainty from the background. These exclusion 
bounds further imply that various popular scenarios appear to be totally incompatible with 
the ATLAS diboson excess. 

One should remark that various possible scenarios giving rise to the observed ATLAS 
excess have been considered so far. While spin-1 resonances have been investigated by 
many authors [8-27], spin-0 and 2 SM-singlet resonances have received far less attention, 
see Refs. [28-31]. Here we go beyond previous studies by setting up the complete effective 
theory for neutral resonances. We also perform a full statistical analysis of the ATLAS 
search, including the extraction of the width from the shape of the excess. 

2 Effective Field Theory for neutral resonances 

In this section we introduce the EFT of SM-singlet resonances of spin 0 (CP even and odd) 
and spin 2 coupled linearly to the SM. We denote the mass of the resonance with m and 
assume that it is much heavier than the electroweak (EW) scale, m 2 v 2 , m^, m?, m\ 
etc, which is an excellent approximation for a hypothetical 2 TeV resonance. 

We will use field redefinitions (or, equivalently, equations of motion) to reduce the 
number of independent operators. The leading interactions will be dimension-5 operators 
and we denote them generically by Ox with coefficients /y 1 where the fx have dimension of 
mass. The region of validity of the EFT is set by the condition that one can neglect higher 
dimensional operators. The most severe restrictions come from operators with additional 
derivatives on cj), such as d 2 (j>G‘j 1L ,, which require us to impose the condition 


( 2 . 1 ) 


m < M. 


where M denotes the cutoff of the theory, at which the nonrenormalizable dimension-5 
operators become resolved by new states of mass M. 


In order to estimate the maximal size of the couplings f x , 
sional Analysis (NDA) which gives 


we can use Naive Dimen- 



( 2 . 2 ) 


Using Eq. (2.1), the maximal allowed size /y 1 is at most of the inverse EW scale for m ~ 2 
TeV. In many UV completions, the coupling is expected to be weaker than the bound (2.2). 
For instance, if the nonrenormalizable coupling 



is resolved in the UV by a heavy fermion of mass M, then one expects 



(2.4) 
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where a s is the strong coupling at the scale M. and the estimate is obtained by taking the 
coupling of cj) to the fermions < 47 t. 

We now list the complete EFT’s for the cases of spin 0 (CP odd and even) and spin 2. 


2.1 Spin-0, CP-even 

The effective Lagrangian for a neutral, CP even, spin-0 resonance reads 

Ah = ^(/g 1 Ku ) 2 + fw 1 ( W P 2 + fs 1 ( B + fv 1 \D^H\ 2 + f- 1 Re(-y t Hi R q L )) , 

( 2 - 5 ) 

where H = IG 2 H . In order to avoid issues with flavor violation, the operators including 
fermions are expected to be roughly proportional to the SM Yukawa couplings, hence here 
we show only the one involving the top quark, denoted by q R and t R . 

A priori, one could have written two more operators (that are also relevant for diboson 
production at the LHC): 

0' H = ^d 2 \H\ 2 , O'lj = (t>\H\ 2 . (2.6) 


The operator O ^ generates a mass mixing after EWSB as well as a tadpole for cj) that 
induces a vacuum expectation value (VEV) for this field (or shifts an existing one). It can 
be eliminated by a held redefinition of </> in favor of O r H , which leaves only kinetic mixing. 
The operator O r H in turn can be eliminated via the Higgs equations of motion in favor of 
Oh and Ot- The resulting Lagrangian (2.5) does neither have mass nor kinetic mixing 
between 4> and the Higgs, nor does it induce any VEV for (j>. We will see an expicit example 
in Sec. 3.1. The operator 0|Lf| 4 gives similar effects as the operator 4>\H\ 2 , but suppressed 
by an additional factor of v 2 /m 2 . 

The decay width resulting from the above Lagrangian (2.5) is then 1 


m 3 f 8 3 1 1 \ 

^\7E + & + J% + m) ' 


(2.7) 


The partial widths can easily be extracted from T, see App. B. 

A brief comment about the operator Ot is in order. The latter can generate couplings 
of cj) to gluons and photons at one-loop. Even though these cannot be written as local 
operators, for our purposes (he., on-shell production) we can represent this diagram by a 
complex contribution to e.g. the 4>GG coupling 2 


Hfc 1 ) ~ /t 1 ^(0.0014 - 0.0044 i ), 


( 2 . 8 ) 


where we have taken m = 2 TeV. Eq. (2.8) can simply be obtained from the corresponding 
expressions of the Higgs couplings to gluons, see e.g. [32, 33]; note the presence of the imag¬ 
inary part due to the tt-mass threshold. Using NDA, 4 tt/t > m, we obtain the estimate 

x We use that mw, mml m 2 . The partial decay width to top quarks is suppressed by a 

relative factor of m 2 /m 2 , see below. 

2 A similar expression can be given for the (j >77 and (jryZ couplings which also receive contributions 

proportional to from the W-loop. 
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| A f G 1 1 < (430 TeV) 1 , and we can safely neglect this contribution to the production of 
cp. Moreover, Ot can induce decays to top quarks with partial width = |which 

is suppressed by a power of m 2 /m 2 compared to the other decay channels. Only if f^ 1 is 
much larger than all the other couplings will this contribution matter. An upper bound 
can again be derived using NDA, one finds that F t j < 70 GeV for to = 2 TeV. 


2.2 Spin-0, CP-odd 


The Lagrangian for a CP odd, spin-0 resonance is 

To- = <t> (/g 1 g%g% + fw 1 w; lu w; lu + /- 1 b iw b, w + /- 1 im( yt Hi R q L )) . (2.9) 

where F )W = \(. pvpa F p(T . The additional operators 

O n = (Pd p i[H^D p H-D p H^H], O t = (f> W), (2-10) 


where ip runs over the chiral SM fermions (ip = u l R , d l R , e l R , l' L and q l L ), can all be eliminated 
by appropriate field redefinitions 3 in favor of Ot- Models giving rise to this effective theory 
have recently been considered in Ref. [29] in the context of the ATLAS results. Notice that 
for Ot, the same comments as in the CP-even case apply. 

The decay width is given by 


TO 


8 


1 


+ t + <2 ' n) 

which is identical to the CP even case, except for the absence of the operator Oh- Notice 
that our results agree with those of Ref. [31] whereas w.r.t. Ref. [29] we find a discrepancy 
of a factor of 4. The partial widht to top quarks is again given by T^ = and can be 

neglected. 


2.3 Spin-2 

We now give the effective Lagragian for CP-even spin-2 fields. A massive spin-2 resonance 
is described by a symmetric-traceless (ST) field <f> pv . As is well known [34, 35], a consistent 
description requires in addition a scalar field (denoted here by x), which enforces transver- 
sality and removes the unphysical longitudinal degrees of freedom, i.e.sets = 0, such 

that only the five physical polarizations remain. Rs equation of motion are algebraic, i.e. x 
is a non-propagating auxiliary field (in the absence of sources it simply vanishes, X = 0). 
We do not write the free Lagrangian here (see however Sec. 3.2) but rather directly give 
the propagator, which in the basis {(p P v,x) reads 4 

P = i (pV n !£ \ 

\d? k i, k -\ s?C c2 + 2m2 '>J 

^Without loss of generality we have assumed the fermion operators to be flavour-diagonal, though 

not necessarily flavor-universal (e.g. ft R ^ fu R )- It should be kept in mind that the degree of flavor- 

nonuniversality is highly constrained by data. In any case we take only the top-Yukawa coupling to be 
nonzero. In making field redefinitions of the chiral fermions, one should keep track of anomalies which will 
generate contributions also to the coefficients of Og,w,b- 

4 Typically the propagator of the massive spin-2 case is given for the reducible representation +ri )ll ,x, 

see e.g. [36]. Here we prefer to display explicitly the decomposition into the irreducible components. 
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with 


(2.13) 


nz = 5^5 V \ - AtSWk 




{p cr} m 2 {p 


j' 3m 4 


biPb^kf h 

IX/ IX/ A/J n 


} 


where the curly brackets denote ST, i.e.Xs^x = \X RV + \X vtl — ^r) lw Xp. In particular, P 
mixes the scalar and tensor degrees of freedom, but only the tensor degrees of freedom have 
physical poles. Notice that on-shell II is simply the projector on transverse, symmetric, 
traceless fields, in particular at k 2 = m? one has II p&k p = 0 etc. As usual, the projector 
can be written in terms of polarization tensors, wich for completeness we collect in App. A. 
The following three observations will further simplify the analysis. 


1. As we are only interested in amplitudes for processes near k 2 = m 2 , only the tensor- 
tensor part of the propagator will matter. In particular, the source of the field x ~ 
which is non-zero in general - will not contribute. 

2. As the tensor-tensor propagator above is transverse on-shell, any source that is just 
a total derivative of the kind dpJ u , (e.g., dpd u \H 2 \) will not contribute near the pole. 

3. Any source that is conserved (such as F^F v y p ) will only receive contributions from 
the term proportional to the identity 

The most general bosonic effective Lagrangian linear in a ST held cf)^ up to dim-5 is then 
simply 


/’bos 

L '2+ 


pv \ JG 


ft 


_ 1 G“ G“ 


pp^up 


+ f - 1 w i W i 

r J W [ip vp 


+ fs l B 


ppB up + DpH^D v H 


(2.14) 


while the most general fermionic source Lagrangian reads 

£ 2 + = 4>pv Im(^D^)j , 


(2.15) 


where the sum over chiral fermions (if) = u l R , d l R , e l R , i l L and q l L ) is understood. We 
remark that unlike in the scalar cases, even the light SM fermions need to be kept, as one 
cannot eliminate them via their equations of motion. Without loss of generality we have 
diagonalized the operators O^, but in principle allow non-universal couplings 
One should keep in mind though that flavor-nonuniversality (e.g. fd R / f SR ) is highly 
constrained by data. It is crucial that one uses the above Fierz-Pauli propagator for the 
computation of the scattering amplitudes arising from Eqns. (2.14) and (2.15). 5 

For the decay width resulting from the above Lagrangian we find (we review in App. A 
the relevant polarization tensors) 


F = ^(A + A + J_ + _j_+M 
^\fh fw fl 12 /V 44 J 


(2.16) 


5 Another consistent possibility would be to introduce Goldstone fields <j> and to render the Lagrangian 
gauge-invariant under linearized general coordinate transformations, and then adopt a gauge in which the 
propagators simplify. However, we stress that this procedure also fixes the sources for the auxiliary and 
Goldstone fields, which cannot be ignored in this case. We will make some more comments on this in 
Sec. 3.2 in the context of a specific example. 
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where denotes the gauge-multiplicity of i/j (QCD plus EW), e.g. N = 6 for a LH 
quark doublet. Our results agree with Refs. [31, 37-39]. We refer again to App. B for the 
decomposition of T into partial widths. 

For completeness, we mention that all symmetric-traceless CP-odd sources up to di¬ 
mension four made from the SM fields are total derivatives of the kind mentioned in point 2 
above, and as such do not contribute to resonant production from dimension-5 operators. 6 
On the other hand, dimension-7 operators are always suppressed by additional powers of 
rn'z W /M 2 , giving only very small cross sections and widths. We therefore do not include 
a CP-odd spin-2 particle in our analysis. 

3 Scenarios 

The purpose of this section is to give a few well-motivated scenarios for the effective theories 
described in Sec. 2. 

3.1 Higgs portal 

Consider a neutral scalar field, of mass m, interacting through the Higgs via the interaction 

C = —^4>\H\ 2 . (3.1) 

The parameter /r has dimension of mass and might itself be an effective interaction resulting 
from some renormalizable coupling g& 2 \H\ 2 after $ obtains a vacuum expectation value 
4> = u + 4>. 

We now make the shift cf) —> </> — \H\' 2 . After this field redefinition the Lagrangian 

becomes 

c = -^d^d^Hl 2 + l -^(d^\H\ 2 ) 2 + ... (3.2) 

mi 2 m- 

<t> 

where the ellipsis denotes unobservable modifications of the Higgs potential. The first term 
is the effective interaction discussed below Eq. (2.5) with the identification 

f'n = ~ ■ (3-3) 

The second term leads to modifications of the Higgs couplings. In order to avoid too-large 
deviations inconsistent with experiment, we will impose that \xv <C m^, which implies 
that v <C f' H . To arrive at our standard basis, we use the Higgs equations of motion (or, 
equivalently, make the field redefintion H —» H + -^ 2 (f)H) to find 

^d 2 \H\ 2 = ^cj>(\D^H\ 2 -y t Re(Hi R q L )+m 2 H \H\ 2 -2\\H\ 4 )) . (3.4) 

m z m z V / 

The last two terms (that results from the Higgs potential V = —m 2 H \H\ 2 + A|i7| 4 ) can 
be neglected, as they are suppressed w.r.t. to the original interaction (3.1) by a factor of 

5 We do not agree with some of the results of Ref. [31]. Their fermionic source = Ini is 

CP even. In fact the CP-odd sources Re< 27 5 7 ( M cb)q and Keq'y( ll ,d v )q are total derivatives. 
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m 2 H /m 2 , \H\ 2 /m 2 * . We then have 


fn = It 


with the remaining vanishing. 


9 

m 

2/x 


(3.5) 


3.2 Spin-2 Lagrangians from warped extra dimensions. 

In this section we derive the massive interacting spin-2 Lagrangian from a warped extra 
dimension [1]. According to the general discussion in Sec. 2, we expect the presence of an 
auxiliary field. Moreover, in the extra dimensional construction we arrive naturally at a 
theory containing Goldstone modes as extra-dimensional components of the metric, which 
one can simply set to zero in a ” unitary” gauge. 

We then consider a 5d compactification in the metric background 

ds 2 = ( kz)~ 2 (dx 2 — dz 2 ) = 7 Mivdx M dx N , (3.6) 

where z denotes the 5th coordinate zo < z < z± and k = z^ 1 the Anti-de-Sitter curvature . 

After decomposition in Kaluza Klein (KK) modes, the kinetic Lagrangian of the fluc¬ 
tuations of the 5d metric becomes' 


£kin = 4>lv{d 2 + m: 




+ \r r (& 2 +mDr,+\x n (s 2 +mDx" - +mDr 


- 18,^ + \e„X n + + m„ x ” + 


(3.7) 


Notice that the Lagrangian is completely diagonal in the KK modes. Here, </>” and cj) n 
denote the Goldstone modes originating from the extra-dimensional components of the 5d 
metric, 8 and is the above mentioned auxiliary field. Setting to zero the Goldstone fields 


cP n = 0, 0” = O, n± 0, (3.8) 

one arrives at the unitary gauge, which is precisely the Fierz-Pauli Lagrangian [34, 35] 
leading to the propagator (2.12). Instead, one could adopt the Feynman gauge, in which 
the terms in the second line of Eq. (3.7) are cancelled by an appropriate Fadeev-Popov 
procedure. In Feynman gauge, the propagators are especially simple, in particular, all 
fields have the same mass and do not mix; observe that the held x becomes propagating 
and has a ’’wrong sign” kinetic term. We however stress that in this case, the sources for 
all fields, 0”^, \ n , 4^ an d 4> n have to be taken into account. In the following we will employ 
the unitary gauge (3.8) in which case we only need to consider the source for the ST held 
6 n 

'We refer the reader to Ref. [40] for details, in particular the precise relation of the various fields 

to the 5d metric. Eq. (3.7) is obtained from Eq. (3.5) of [40] by use of the 5d wave functions f s = 

x/2z s Js (m n z)/zi J 2 (m,nZi), where the J„ denote Bessel functions. The masses are solutions to Ji(rri n z\) = 0 

and x have wave functions /o, has wave function /i, and cj> has wave function /o. 

8 The field (/> M does not have a zero mode, while (j> has a zero mode that is not eaten and corresponds to 
the radion. 
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As for the interactions, typically two scenarios are considered. In the brane model, 
all SM fields are localized on the IR brane [1], while in the bulk model they are allowed 
to propagate in the bulk [38, 41, 42], In the latter case, the gauge fields have flat 5d 
profiles, the RH top and the Higgs fields have profiles peaked towards the IR brane, and 
the remaining matter fields have profiles that are flat or peaked towards the UV brane. 9 For 
our purposes it is good enough to approximate the bulk model by IR brane localized RH 
top and Higgs fields and completely ignore the other quarks and leptons. The interaction 
terms for IR-brane localized fields are given by 

£ matte r = ( 2 _ Im i R ^D u t R ) , (3.9) 

where Mp is the reduced Planck mass. In the scenario with all SM fields localized on the 
IR brane, there are identical contributions to the remaining SM fermions. Gauge fields 
couple as (for any n/0): 


^gauge _ z lk jjn 1 7 zr /- _ r l T 2 X n [1 + J2{x n ) ] 

int ~~ Mp U > u - ro + ri + v 


(3.10) 


where x n = z\m n and the quantities ro and r\ denote possible brane kinetic terms (BKT) 
[43], and V = log(A;zi) ~ 36 is the volume of the extra dimension. An IR brane-localized 
gauge field is described by the limit r 1 —> oo, or £ n = l. 10 

For the bulk model, the effective Lagrangian for the first KK mode of mass m is then 
given by £ 2 + defined in Eqns. (2.14) and (2.15), with the couplings 

= f G = f w = f B = --( cm)- 1 , ft R = -- x r 1 (3.ii) 

Zitv tv tv 


where x\ = 3.83, and k = k/Mp the RS coupling parameter. For the brane model, one has 
instead 

fH = y~_Xi 1 ’ fa = fw = Sb = -—ah 1 , U = -—xJ 1 . (3.12) 

Zitv tv tv 

According to our general formula Eq. (2.16), the bulk and brane model’s total widths are 
respectively 


Fblllk — 


(13 + 144£i) (ftxi) 2 m 
9607T 


Fbrane — 


283(kxi ) 2 m 
9607T 


(3.13) 


In the RS bulk model, the terms proportional to contribute 0.2%, 26%, and 56% to the 
total width for r\/V = 0, 0.2 and 0.5 respectively. 

9 We remark that such a scenario features other states, typically lighter than the KK graviton, whose 
phenomenology will severley constrain the model. We will not further consider these model-dependent 
constraints in this work. 

10 We remark that in this limit the KK modes of the gauge fields become strongly coupled (/k K ~ g 2 (V+ri), 
hence to avoid the non-perturbative regime one would demand \/V + n < 47r/g. Note also that the gauge 
KK modes decouple from the IR brane in this limit. 
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3.3 Radion/Dilaton 

In the warped extra-dimensional scenario considered in the previous section, the field />° 
corresponds to the radion which describes the fluctuations of the size of the extra dimen¬ 
sion. It is massless in the background (3.6) but by a suitable stabilization mechanism it 
acquires a mass [44-46]. Athough its five-dimensional wave-function is deformed by the 
stabilization mechanism 11 we will assume that these effects are small and its couplings are 
thus approximated by those of the massless case. One finds (see e.g. Ref. [49]) 

Ant = Gco + ro + n) + 2 |D/tjf/|2 ~ 2Re ^ ytHtRq ^) (3 ' 14) 

There is an additional coupling proportional to the Higgs potential V(H) = — m 2 H \H\ 2 + 
A|i7| 4 . Eliminating the operator 4>\H\ 2 will result in negligible corrections to the operator 
coefficients Z# 1 ,// 1 suppressed as m 2 H /m 2 . 

It is customary to treat the radion interaction scale defined as 


/rad 


\/6 Mp 
kz\ 


as a free paramter. In the bulk model one then finds the couplings 


fH = fr 



/g = fw = Ib = 4V / ra d 


(3.15) 


(3.16) 


where we have assumed r* <C V. The brane model is again obtained by sending r\ —> oo. 
The couplings to the gauge boson field strength vanishes in this case, and one is left with 
only 


fn = It 



(3.17) 


Interestingly, the brane model effective Lagrangian precisely conincides with the Higgs 
portal scenario with the identification / ra d = m 2 //r. In either case, the field <f> just inherits 
the Higgs couplings suppressed by a factor u// ra d- As the couplings to gauge bosons are 
always small, the decay width comes entirely from fn in both models 


T = 


m 


8 irf 2 


rad 


(3.18) 


As explained in Sec. 2.5, the decay to tops are suppressed by m 2 /m 2 and do not contribute 
to the total width. Finally we recall that the radion is closely related to the dilaton of 
nearly conformal extensions of the SM, so that very similar results hold in this case. 


4 Characterisation of the ATLAS diboson excess 
4.1 Data, background and local significances 

The ATLAS collaboration has recently presented a search for narrow resonances decaying 
to electroweak bosons with hadronic final states using the 8 TeV LHC dataset [7]. This 

11 See Refs. [47, 48] for some analytic expressions for the stabilized profile. 
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fl r 

b r 

h 

b 

WW 

4 

7 

2 

9 «7+0- 42 
Z - D ‘ -0.40 

1-84+1? 

1 30+ 0 - 28 
l.OU —o 24 

13 

6 61 +1,22 
D - Di -1.01 

WZ 

5 

8 

2 

q I o+0.50 
°* lz -0.42 

2-08+1 

i a 1+0.33 
i -^ i -0.26 

15 

c D-I+1.04 
^■81—0.95 

ZZ 

5 

3 

1 

n qi+0-23 
u ' yl -0.20 

0-55+1® 

n , ?4+ 0 - 12 
u - o ^-0.09 

9 

-I O+0.50 
1 -°-0.42 


Table 1. Data and background obtained from [7] in the three bins r = {[1.85,1.95], [1.95,2.05], 
[2.05, 2.15] TeV} and in the whole excess region of the dijet mass spectrum, for the WW, WZ and 
ZZ selections. 


dataset has 20.3 fb^ 1 integrated luminosity. The weak bosons from massive resonances are 
highly boosted and are thus reconstructed as a single jet with large radius using advanced 
reclustering, grooming and filtering algorithms. The expected background is dominated 
by dijets events from the QCD background, which is huge but does not feature potential 
resonance structures. 

Boson-tagging cuts are applied to the selected dijet events, asking for subjet momentum- 
balance and low number of associate charged particles tracks. Each jet is then tagged using 
a narrow window on the jet mass rrij, asking for rrij to be close to the W or Z mass. In 
the analysis, a jet is identified as a IE if rrij £ [69.4, 95.4] GeV, and is identified as Z if 
rrij £ [79.8,105.8] GeV. The W and Z masses being close, these two ranges overlap. There 
are thus three disjoint tagging regions, that we label as IE-only, IE or Z (noted W/Z), and 
Z-only. 

A local excess of observed events appears in the dijet spectrum near 2 TeV. The num¬ 
bers reported in Ref. [7] (and its extra material [50]) in the three bins rrijj £ [1850,1950], 
[1950,2050], [2050,2150], that we refer to as the excess region, are shown in Tab. 1. The 
expected dijet background in each bin is also shown. The background is partly determined 
from a fit to the whole dijet spectrum, and is thus subject to some uncertainty. 

As a first step, one should check the statistical significance of this excess. Assuming 
Poisson statistics for the observed events in each bin, we first compute the p-value of a 
discovery test in every bin. This computation is done with and without taking into account 
the background uncertainties, that we model using a nuisance parameter 6 £ [0 a ,0b] with 
a flat “prior” distribution. 

The likelihood for one of the bins r simply reads 

Q\ _ ( S r + K + 0) h r e -Sr-br-0 
S r ,U) = 7 j 

n r \ 

The nuisance parameter is eliminated by maximising this likelihood with respect to 8 for a 
given s r , L(s r ) = maxg L(s r , 0). The statistical significance Zq for the existence of an excess 
is obtained by computing the probability density f q for q = — 2 log[Z(s r )/ max Sr Z(s r )] and 
evaluating the observed p-value p = f°° dqj Q . The p-value is further translated into a 
standard significance by Zq = <J>~ 1 (1 — p), where is the standard cumulative Gaussian 
distribution. One allows for both upward and downward fluctuations. The significance of 
this discovery test is computed for each bin. The values, shown in Tab. 2, typically go 
beyond two sigmas in the central bin. 


(4.1) 
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Zo 

without syst.) 

^0 

(with syst 

) 

B 0 

WW 

< 1 

2.8 

< 1 

< 1 

2.5 

< 

1 

0.9 

13 

0.8 

WZ 

< 1 

3.0 

< 1 

< 1 

2.7 

< 

1 

1.0 

22 

0.8 

ZZ 

2.8 

2.1 

< 1 

2.6 

1.8 

< 

1 

21 

6 

1.1 


Table 2. Local discovery tests for the [1.85,1.95], [1.95,2.05], [2.05,2.15] bins of the dijet mass 
distribution with WW, WZ and ZZ tagging. Left and middle columns: discovery significance 
without and with systematic uncertainties. Right column: discovery Bayes factor. 


We also introduce a Bayesian discovery test, defined as 

p{h r \s r ) 


Bn = 


p(n r \s r = 0) 


This expression takes the simple form 


£n = 


f ds r L(s r )Tr(s r ) 

W) 


(4.2) 


(4.3) 


It turns out that the prior for the signal ir(s r ) is entirely fixed from general considerations. 
Indeed, the measurement being a counting experiment, we already know a priori that b r +s r 
follows a Poisson distribution. The parameter of this Poisson distribution has to be chosen 
to be b r , which is known a priori, in order not to bias the discovery test. This then fixes 
7r(s r ) to be 12 


Tr(s r ) 



(b r + s r )! 


(4.4) 


The values of the discovery Bayes factor are shown in Tab. 2. One can see that the values 
of Bq are beyond the threshold of moderate evidence for the central bin. 

It follows that both frequentist and Bayesian discovery tests provide a moderate evi¬ 
dence for the existence of a local excess over the QCD dijet background. We conclude that 
this excess is significant enough to deserve attention, so that we proceed in the analysis. 


4.2 Mass and width reconstruction 

As the data are provided in several bins, it is possible to analyse the shape of the hypothet¬ 
ical signal. Even though the statistics of the excess is fairly low, we emphasize that there 
is no reason that prevents to apply a rigorous shape analysis. Whether or not the data are 
informative enough should be decided by the outcome of the analysis. Notice that, as the 
excess is observed in more than one bin, one can expect both an upper and lower limit on 
the width of the resonance. 

In what follows the bins of the rrijj distribution are labelled by the index r. Contrary 
to the analysis on the total event numbers, here we do not combine the events of the 

12 In [51] it will be shown that this particular prior provides a good connexion between discovery Bayes 
factor and frequentist statistical significance. Also, notice that we do not implement the background 
systematic error in the Bayes factor. This is because this type of systematic uncertainty approximately 
cancels out in the Bayes factor, as will be shown in [51]. 
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three selections WW, WZ, ZZ, and rather perform the shape analysis for each selection 
separately. It will be clear from next section that a more evolved analysis combining the 
three selections would bring only little extra information. 

The likelihood containing the shape information appears naturally from the full likeli¬ 
hood L = L r , by factoring out the likelihood for the total event number, L = £ tot £ sha P e . 
Explicitly, the shape likelihood reads 


^shape 


n 

r (bins) 



(4.5) 


Note that the factorisation L = / y tot L shape makes clear that a shape analysis of the diboson 
excess is truely complementary from the total event number analysis, because each analysis 
rely on mutually exclusive pieces of information. 

We denote the shape of the expected signal by a distribution f m .. normalised to one 
(i.e. a density). The shape of the signal is modelled assuming a resonant amplitude, and 
the background is assumed to be flat near the peak of the resonance. The narrow-width 
approximation is assumed, i.e. one takes T/m -C 1, that will be well verified a posteriori. 
Given these standard assumptions, the rrijj distribution is then distributed following a 
Breit-Weigner shape, 

f m n (77)2 j _ m 2y _|_ m 2p2 ( ) 

The expected content of the bins is obtained by integrating over this distribution, 


n r — ntot I f m jj ’ (4-7) 

J bin r 

an one will note f r the shape density integrated over a bin, f r = f hin r dnijj fmjj ■ We 
consider the three bins centered around 2 TeV, and assume no signal event elsewhere. 

We also take into account the systematic uncertainties relevant for the shape of the 
signal. These are the uncertaintites on the jet reconstruction (see [7]), that tend to smear 
the resonance shape. The sources of error are the jet px resolution, the jet px scale and 
the jet mass determination, associated respectively to the nuisance parameters d res , d sca i e , 
<5 m , affecting the rrijj mass. The magnitude of these errors is small with respect to one, so 
that they can be written in the linear form 


( 1 -j- (5 r es T d'scale T dm) • (4.8) 

All these uncertainties are modelled using Gaussian nuisance parameters d res , dscaie; d m 
with zero mean and respective standard deviation <r res = 0.033, (T SC ale = 0.02, cr m = 0.03 
(see [7], Tab. 4). 

These three nuisance parameters being independent, they can be rigorously combined 
into a single Gaussian nuisance parameter 6 with zero mean and variance given by 

^ = ^res + Gscale + a m ■ ( 4 - 9 ) 
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Mass 

68 % CL 

[GeV] 

95% CL 

Width 

68 % CL 

[GeV] 

95% CL 

WW 

[1895,2091] 

[1797, 2190] 

[39, 99] 

[22,150] 

WZ 

[1895,2091] 

[1797, 2189] 

[42, 98] 

[26,144] 

ZZ 

[1852, 2047] 

[1755,2145] 

[15,71] 

[6,139] 


Table 3. One-dimensional confidence intervals at 68% and 95% confidence level for the mass (to) 
and the width (r) of the hypothesized resonance. The intervals are computed independently for 
each subchannel. 



m( GeV) 


m( GeV) 


m( GeV) 


Figure 1 . Probability densities in the m — T plane, for WW, WZ and ZZ selections (from left to 
right). The green, yellow and gray regions correspond respectively to 68%, 95% and 99% confidence 
level. 

The event number in a given bin depends thus on m, T, <5, so that the complete likelihood 
for the shape analysis of the diboson excess reads 

L(m,T) = f dd [f r (m,T,6)f r tt(S) (4.10) 

r (bins) 

For the mass and width of the hypothesized resonance, one assumes log priors vr(m) oc m -1 , 
7r(r) oc r -1 , which are the most objective priors for dimensionful quantities. The confidence 
regions are drawn from the posterior density, which is given by p{m , T) = L(m, r)7r(m)7r(r). 

The one-dimensional confidence intervals for mass and width are given in Tab. 3. The 
systematics errors increase the mass CL bounds by roughly ~ 5% and the width CL bounds 
up to ~ 20%. Using a flat prior instead of a log prior changes the bounds by roughly ~ 10%. 
The two dimensional confidence regions in the m — T plane are shown in Fig. 1. In the 
following, we shall quote the results from the WZ selection, which contains the largest 
event number. 

5 Statistical analysis of the diboson rates 

Having studied the shape of the diboson excess, we now turn to the analysis of the overall 
event numbers, i.e. the total rates over the excess region. The likelihood analysis for a 
set of overlapping selections is a somewhat unusual exercise to carry out, so that we shall 
provide a detailed explanation of the statistics involved. For clarity, in the following we 
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will use rigorous probability notation. The hypothetical event number in a given selection 
is taken as random variable, denoted by N. Specific values of event numbers are denoted 
by n, and P(N = n) is the probability of N for taking the value n. The expected event 
numbers are denoted by A, and the observed event numbers are denoted by n. 

5.1 The statistics of hadronic weak-boson tagging 

The mass distribution of a fat jet coming from a IT or Z is peaked at the boson mass, m\y 
or mz- The jets can be therefore tagged as IT and Z by requiring rrij to be close to mw,z- 
In the analysis of [7], a jet is identified as a IT if rrij E [69.4, 95.4] GeV, and is identified as 
Z if rrij E [79.8,105.8] GeV. The W and Z masses being close, these two ranges overlap. 
This implies there are three disjoint regions to tag the jet: 

• If rrij E [69.4, 79.8], jet is IT-only, 

• if rrij E [79.8,95.4], jet is IT or Z (noted W/Z ), 

• if rrij E [95.4,105.8], jet is Z-only. 

These tagging regions will be labelled by I. The expected rrij distributions have been 
provided in Fig. lc of [7]. These distributions as well as the tagging regions are shown in 
Fig. 2. 

Note that the distributions for true IT and Z have been generated assuming a bulk 
RS KK graviton signal. From Sec. 3.2, it is clear that f^ 1 fy \ so that the bulk RS KK 
graviton decays mostly to longitudinally polarized IT and Z. However, the weak boson 
widths being narrow and the final shape being strongly widened by the detector effects, we 
expect the IT, Z distributions of Fig. 2 to hold for any polarisation of the weak bosons to 
a very good approximation. 

Using the distributions of Fig. 2, it is possible to estimate the tagging probabilities, 
given one of the two hypothesis for the underlying true boosted particle, {True W, True Z} 
that we will label by X. What we compute is thus the conditional probability p(I\X). The 
conditional probabilities for tagging a true IT and a true Z are computed from Fig. 2 and 
shown in Tab. 4. These numbers are consistent with the ones found in [21]. 


p(I\X) 

IT-only 

IT or Z 

Z-only 

True IT 

0.253 

0.366 

0.034 

True Z 

0.112 

0.398 

0.211 

True j 

0.025 

0.035 

0.023 


Table 4. Conditional probabilities p(I\X) for IT, Z and background jet tagging. 

Moreover, fat jets can also arise from the QCD interactions. The distribution for a jet 
coming from the QCD dijet background has been simulated in [7] (see Figure 1 c there), and 
appears to be nearly constant over the tagging regions. Using the simulated distributions, 
we can deduce the probabilities for mis-tagging a jet from the QCD background as a weak- 
boson jet. Finally, the total probability for tagging a IT, Z or j as a weak boson V is just 
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Figure 2. Mass distributions of a jet arising from a W (blue curve), Z (black curve) and the 
QCD dijet background (red curve). Blue and black dotted lines represent the W and Z tagging 
regions respectively, giving rise to the three disjoint tagging regions, IT-only, W/Z and Z- only. 


Sk _^ 


Nio 

Nu 

N 0 i 

W -only 

W/Z 

Z- only 


Figure 3. A picture summarizing the event numbers over the various weak-boson tagging regions. 


obtained by summing the probabilities over the three region. One gets P{V\W) = 65%, 
P{V\Z) = 72%, P(V\j) = 8%. 

Before closing this subsection, it is instructive to focus on the counting statistics for 
the tagging of a single jet. This part can serve as a statistical toy-model for the upcoming 
analysis of the diboson excess. Indeed, most of the ingredients for the diboson analysis are 
already there, though applied to a simpler problem. 

Let us denote the tagging regions VT-only, W/Z , Z-only as 10, 11, 01, and labelled 
by / G 10,11,01. The first number of the region name means that the region potentially 
contains a W if equal to one, and does not contain a W if equal to zero. The second 
number of the name works similarly for the Z. These notations will be convenient later. 

The event numbers in each of these regions are denoted IVio , -/Vn, JVoi. These events 
follow independent Poisson statistics with parameter Aio, An, Aoi, 

\ ni p-\i 

P(N I = n I \X I ) = 1 , ■ (5.1) 

Assuming an expected event number Xx = (Ayy,Az) for the true W and Z, the A/ are 
expressed as 

\i = ^P(I\X)X x . (5.2) 

x 
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Equations (5.1),(5.2) put together provide P(Nj = n/jAx), the probability of observing nj 
events in the region I for given expected event numbers Ax- Taking this probability as a 
function of Ax provides the likelihood function for Ax, for an observed event number nj. 

Let us now assume that only the number of events that contain all possible IT-tags 
and all possible Z -tags are reported. These numbers are defined as 

Nw = Ni o + Nn , Nz = Nqi + Nn . (5.3) 


This configuration is pictured in Fig. 3. Clearly, the statistics of N\y and Nz are not 
independent, because of the common region 11 where the jet is either W or Z. Rather, the 
Nw, Nz follow a bivariate Poisson statistics, given by 


P(Nw = nw,Nz = nz\Xi) 


E 

niO+ n ll= n w 

noi+nn=n z 


\noi \nio \nn 

A 01 A 10 A 11 c -Aiq-Aqi-Aii 
n 0 i!ni 0 !nn! 


(5.4) 


The mean of (Nw, Nz) is given by (Aio + An, Aqi + An), and the covariance matrix is 


Aio + An An 
, An Aqi + An 


(5.5) 


Plugging Eq. (5.2) into Eq. (5.4), one gets the probability of getting (nw,nz) events for 
given expected event numbers Ax- Taking this probability as a function of Ax provides 
the likelihood function for Ax, for an observed event number nw and nz- 


5.2 Statistics for the ATLAS diboson excess 

The probability for the tagging of two fat jets are obtained by combining the probability of 
tagging a single jet, see Tab. 4. For the tagging of two jets, six tagging regions are obtained, 
by combining the labels W-only, W/Z and Z-only in all inequivalent ways possible. The 
index I of the tagging regions takes then the values 

I € {(W, W), (W, Z), (Z, Z), (IT, W/Z), (Z, W/Z), (W/Z, W/Z)} . (5.6) 

The true events can be either a pair of weak bosons, a QCD jet mis-identified as a weak 
boson or two QCD jets mis-identified as weak bosons. The list of the hypothesis of true 
events, is then 

X € {(W, W), (W, Z), (Z, Z), (j,j), ( W,j), ( Z,j )} , (5.7) 

where j stands for background jet. The conditional probabilities P(I\X) are given in Tab. 5. 
The numbers for true WW, WZ,ZZ are consistent with the ones reported in [21]. The 
dijet background corresponds to the true event jj. Pileup effects are assumed to be small, 
so that we do not consider the possibility of having true events as (W. j), (Z,j). On the 
other hand, one may consider a new physics signal giving rise to a IT and a jet or a Z and 
a jet. We include therefore the probabilities P(I\Wj s ), P(I\Zj s ) in our table, assuming 
that the distribution of this signal jet j s is roughly the same as from a QCD jet. This case 
will not be considered in the rest of this work, as the decay of singlet resonances does not 
give rise to such signal. 
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P(I\X) 

(IF, W) 

(W,Z) 

(Z,Z) 

(IF, W/Z) 

(■ z,w/z) 

{W/Z, W/Z ) 

True WW 

0.064 

0.017 

0.001 

0.185 

0.025 

0.134 

True WZ 

0.028 

0.057 

0.007 

0.142 

0.091 

0.146 

True ZZ 

0.013 

0.047 

0.045 

0.089 

0.168 

0.158 

True jj 

6.25 ■ 10“ 4 

11.5 ■ 10“ 4 

5.29 ■ 10“ 4 

17.5 ■ 10 -4 

16.1 ■ 10“ 4 

12.3 ■ 10“ 4 

True Wjs 

6.33 ■ 10“ 3 

6.7- 10“ 3 

0.78- 10“ 3 

18.0 ■ 10“ 3 

9.61 ■ 10“ 3 

12.8 ■ 10“ 3 

True Zj s 

2.80 ■ 10“ 3 

7.85 ■ 10“ 3 

4.85 ■ 10“ 3 

13.9 ■ 10 -3 

16.5 ■ 10“ 3 

13.9- 10“ 3 


Table 5. Conditional probabilities p(I\X) for tagging WW, WZ, ZZ true events and a pair of 
background jets jj. The probabilities for tagging Wj s and Zj s are also included. 


The number of events Nj in each of the disjoint tagging regions / follows a Poisson 
distribution with parameter A /, which is related to the expected number of true events 
(i.e. events before tagging) as 

X I = Y / P(I\X)Xx. (5.8) 

x 

The background expected event number A jj will be obtained later on from the ATLAS 
analysis, once we know the statistics for the events. The A ww, ^wz, A zz are assumed 
to come only from the signal, i.e. the SM diboson background is neglected, following the 
ATLAS analysis. The WW, WZ, ZZ expected event numbers are related to the total 
cross-sections by 

Ax = exB x £ crx , (5.9) 

where C = 20.3 fb^ 1 is the integrated luminosity of the 2012 run, B\yw = B^, Bwz = 
B\yBz, Bzz = B 2 Z where Bw = 67.6%, Bz = 69.9% are the hadronic branching ratio of 
the weak bosons. 

The efficiencies e' x for selecting and tagging the signal are reported in [7], Fig. 2b. One 
gets roughly e' x ~ 0.10,0.13,0.09 with about 20% of relative uncertainty. Note that these 
efficiencies are obtained assuming particular models. 13 Slightly different efficiencies can 
be expected for different spins and couplings. This model-dependence should be taken as 
an extra systematic uncertainty on the efficiencies. As the weak-boson tagging probability 
based on the jet mass is already taken into account through the P{I\X), it has to be 
removed from the e' x by dividing by P(V\X). The efficiencies ex we will use are therefore 
given by 

_ e 'ww _ e 'wz _ e 'zz /c i n'i 

wn P{V\W) 2 ’ P(V\W)P(V\Z) ZZ P{V\Z) 2 ’ 1 j 

so that ex ~ {23%, 28%, 17%}. 

In the ATLAS note [7], the expected event numbers A/ on the disjoint tagging regions 
are not reported. Rather, only the number of events that contain all possible kFIT-tags, 
WZ -tags and ZZ -tags are quoted. We denote them by N\yw , ^wz, Nzz■ It is convenient 

ll Tor example, the bulk RS graviton used for the spin-2 simulation and treated in Sec. 3.2 features the 
couplings /jj 1 3> fy 1 , so that it decays mostly to longitudinal polarisations. 


18 

















AW 

Non 

Nooi 

z,w 

z,w/z 

z,z 
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Nm 

Non 

w,w/z 

w/z , w/z 

z,w/z 

Nioo 

Nuo 

Noio 

w,w 

z,w/z 

z,w 


Figure 4. A picture summarizing the event numbers over the various diboson tagging regions. 

to label the tagging regions with respect to their contribution to one or several of these 
reported rates. 

The labels are shown in Fig. 3. Using this parameterisation for the events, the observed 
events read 


Nww = N W o + iVno + Nm i 

= Nqio + Nuo + iVon + Niu , (5-11) 

Nww = iVooi + iVon + N in ; 


Clearly, these events are not independent. They rather follow a trivariate Poisson statis¬ 
tics [52], 


P(Nww = nwwi Nwz = nwz, Nzz 


nzz\M) = ^TT 
v i 


nj i 


A i'- 1 e 


nil 


(5.12) 


where the nj = {nioo, noio,nooi, nuo, non, nni\ are positive integers running over the 
domain 


V = V{nwwi n wz, nzz) = < 

/ 

ni 

nioo + nuo + nm = nww ) 

noio + «no + non + nm = nwz i 


< 

nooi + non + n m = n Z z ) 


(5.13) 


The mean of this distribution is given by 


Nww = Aioo + Ano + Am 
Nwz = Aoio + Ano + Aon + Am 
Nzz = Aqio + Aon + Am 


The covariance matrix is given by 

(Nww Ano + Am Am 

Nwz A 0 n + Am 


V 


Nzz 


(5.14) 


(5.15) 
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The likelihood associated with the measured values of hww , ftwz, ftzz is obtained by 
taking Eq. (5.12) as a function of the hypothesis ( i.e. A/) and using Eqs. (5.8). Dropping 
an irrelevant constant factor, the likelihood is a function of the the various event numbers 
before tagging Ax (recall that X = {WW,WZ, ZZ,jj}), 

L ( Ax) = En (E P ( J I X ) A *) ; (5. 16) 

v 1 x 

where the observed event numbers appear through the domain V = T>(nwwiftwz,ftzz) 
and nowhere else. The A ww, wz, zz from the new physics signal are further related to the 
total production cross-sections by Eq. (5.9). The evaluation of the expected event number 
from dijet background A jj is discussed in the next subsection. 


5.3 Consistency checks and the background likelihood 


As a consistency check of our analysis, we can verify whether the event numbers tagged 
over the full range m,jj E [1,3.5] TeV in the observed sample are consistent with the 
tagging rates determined in Tab. 5 and with our statistical model leading to Eq. (5.14). 
The observed number of events for the WW, WZ, ZZ selections are given in Tab. 8 of 
[7]. These are fiww = 425, fiwz = 604, hy) l z = 333. Regarding the expected rates, the 
complete region being overwhelmed by the dijet background, we can neglect the signal to a 
good approximation, so that the contributions to all tagging regions are simply proportional 
to the dijet expected event number over the full range, A 1 )] 11 . The ratios of the expected 
event numbers in the WW, WZ, ZZ selections are then obtained using the tagging rates 
Tab. 5, Eq. (5.8) and Eq. (5.14). It comes nww/ n w l Z = 0-63, n i Z z/ n wz = 0-57 which are 
in agreement with the observed ratios within ~ 10%. 

The statistical error on the ratios of the n fu11 being roughly about 10%, this consistency 
check seems to be fulfilled within one standard deviation. However, this naive observation 
is too optimistic, because the event numbers of the WW, WZ, ZZ selections are actually 
strongly correlated. The joint statistics of the three selections is a trivariate Poisson, 
already described above, that now describes the whole dataset (i.e. rrijj E [1,3.5] TeV). 
The mean and covariance matrix are thus given as in Eqs. (5.14), (5.15). The covariance 
matrix reads 


T/full _ \full 

V b ~ A jj 


/3.60 2.98 1.23\ 
2.98 5.74 2.84 , 

\1.23 2.84 3.37/ 


(5.17) 


where one used the values of P(I\jj) obtained in Tab. 5. The event numbers n fu11 being 
large, one can adopt the Gaussian approximation so that the likelihood reads 


Lt\ AJT) = exp 


\ E(^ uU - P(i\n)^)(vt ll )l}(nT - r(An)^f) 


(5.18) 


The maximum likelihood gives —2 logT(A^ 1 ) 11 ) = 10.6. This value can readily be interpreted 
as a compatibility test, whose statistics is a chi-square distribution with 3—1 degrees of 
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freedom. The equivalent statistical significance obtained is Z = 2.6. The compatibility 
is thus lower than the 1 <t deviation naively found when neglecting correlations. This 
level of compatibility can nevertheless be considered as acceptable for high-energy physics 
standards, so that we pursue our analysis. 

After these preliminary sanity checks, we now aim at building a consistent likelihood 
for the dijet background event number A jj over the excess region [1.85, 2.15] TeV. The shape 
of the dijet background has been been estimated in [7] using a smoothly falling distribution 
fitted to the observed dataset over the rn 33 £ [1,3.5] TeV range. A different fit is done 
for each of the three selections WW, WZ, ZZ. To the best of our understanding, each of 
these fits should give close results, because the only difference between the selections lies 
in the m 3 ranges selected. Comparing the m 3 intervals with the slope of the m 33 shape, 
it appears that only a slight decrease with m 3 of the efficiency of the boson-tagging cuts 
might be expected when going from the VC-only to the Z-only region. 

The outcome of the three fits can be seen in Tab. 6 and in Fig. 5 of [7]. Comparing the 
central values obtained from the various fits using the quoted error bars, it appears that 
these fits are compatible with each other only within roughly three standard deviations. 
Again, this naive comparison does not take into account the correlations, i.e. it assumes 
that the fits are independent from each other. These fits being partly based on the same 
dataset, their outcome are actually correlated, which implies that the actual uncertainty 
is smaller than the one naively expected. This implies that the compatibility between the 
fits is worse than what naively expected. 11 The shape systematic uncertainties evaluated 
in [7] are found to be small, so that they cannot help solving this discrepancy. 

In order to establish the dijet background likelihood using both a consistent and con¬ 
servative approach, we shall (i) take the correlations among the fits into account and (ii) 
assume somewhat larger uncertainties than the ones quoted in [7]. The likelihood for the 
expected dijet event number in the excess region [1.85, 2.15] TeV before tagging is approx¬ 
imately given by 15 


Tfe(Ajy) = exp 


\ - p (I\mn)(Vb)u(bj - P(j\jj)\ 33 ) 


(5.19) 


where the P(I\jj) are given in Tab. 5 and the expected values of bj obtained from the fits 
are given in Tab. 1. The covariance matrix 14 is proportional to Eq. (5.15), 


V b 


a ^jj 


/3.60 2.98 1.23' 
2.98 5.74 2.84 
\1.23 2.84 3.37 


• 1CT 3 


(5.20) 


where a is a parameter that we tune to obtain a reasonable level of compatibility between 
the fits. As a criterion for the compatibility, we ask that — 2 log L b (Xjj) be equal to the 

14 We point out that in the current version of [7], the background model formula together with the 
best-fit values for the shape parameters do not reproduce at all the background curve shown in the plots. 
This presumed inconsistency is another motivation to adopt a conservative approach for the background 
likelihood. 

15 This is obtained by neglecting the shape systematic uncertainties and assuming standard error propa¬ 
gation. 
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level of compatibility obtained between the selections over [1,3.5] TeV (see above). The 
criterion is thus 

— 2logLb(Xjj) ~ 10.6 . (5.21) 

The maximum of the background likelihood L^Xjj) is found to be 

A jj = 1422 . (5.22) 

Besides, the coefficient satisfying the criterion Eq. (5.21) is 

a « 0.377. (5.23) 

This fixes the overall amount of uncertainty, so that the background likelihood is fully 
determined. The confidence intervals for A jj are found to be 

A jj € [1151,1756] at 68%CL , [935, 2163] at 95%CL . (5.24) 

Note that the 68 % range translates as error bars Ig'gg, ^' 55 , 91 011 the expected back¬ 

ground events bww , bwz and bzz respectively. As expected, these errors are larger than 
the ones quoted [7], that are shown in Tab. 1. 

In order to model the systematic uncertainty on the background, the likelihood will 
be included into the full likelihood, Eq. (5.28). 

5.4 Analysis of total rates 

In the previous subsections, we have gradually derived the total likelihood that should be 
used to analyse the ATLAS diboson excess. It is given by the product of the likelihood 
derived from the counting statistics, Eq. (5.16), times the likelihood constraining the back¬ 
ground, given in Eq. (5.19). In addition, as noted in [21], information on the counting of 
WW + ZZ and WW + WZ + ZZ are available in the additional material of [7]. The values 
hww+zz = 17, hww+wz+zz = 17 are reported. This introduces two new constraints on 
the event numbers nj of the disjoint tagging regions, 


moo + moi + mio + mil + mil = hww+zz 
moo + wool + mio + mio + mil + mu = hww+wz+zz, 


(5.25) 

(5.26) 


that have to be added to the previous constraints already contained in V. see Eq. (5.13). 
It turns out that only three combinations are allowed, so that the domain V is given by 
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(5.27) 


These numbers agree with the ones reported in version 3 of [21]. The final likelihood 
that we shall use to constrain the cross-sections for a hypothetical signal ay, with Y = 
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Figure 5. Posterior densities for the WW, ZZ expected event numbers before tagging. The 
systematic uncertainty from the background determination is included in the right pannel. The 
Xwz expected event number is set to zero. The 68%, 95% C.L. regions are shown respectively in 
green and yellow. 


{WW, WZ, ZZ} are thus 
L(uy ) — j d\jjLi,(Xjj) 

Y n ( P(I\Y)e Y B Y <TY + P(J\Y), Y B Y a Y +P(J\n)X n 

v> I ' Y ' 

(5.28) 


The expected event number for the dijet background Xjj is constrained by Lb(Xjj) and 
is treated as a nuisance parameter. The term in the second row agrees exactly with the 
likelihood used in version 3 of [21], 

Our interest being in neutral resonances, one should first compare the H{Xwz = 0) = 
{Xww % 0, A zz 0, Xwz = 0} hypothesis with the general hypothesis H = {Xww / 
0, A zz / 0, Xwz / 0}. A consistent way to carry out such hypothesis testing is to compute 
the Bayes factor 


B(Xwz = 0) 


p(H(Xwz = 0)|data)) 
p(i^|data) 


/ dXww,zz 7r(Ay)L(Ay) 1 
/ dXy 7r(Ay)L(Ay) 7 t(A wz) 


Vz=0 


(5.29) 

For the prior of the Ay, as described in Eq. (4.4), we use Poisson distributions with the 
Poisson parameter identified as the expected number of background events by, i-e- P(by + 
Ay |Ay). These priors arise from physical considerations and are conservative as they favour 
the background-only hypothesis. We find 


B(^wz = 0) = 0.96 , 


(5.30) 


which implies that the A wz = 0 hypothesis is essentially as credible as the A wz % 0 
hypothesis. This conclusion remains true for the A wz = 0 and A zz = 0 hypothesis as well. 
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Figure 6. Constraints from the ATLAS diboson excess [7] on a generic neutral resonance with 
spin 0. A total width of 150 GeV and 20 GeV is assumed on the left and right pannel respectively. 
Purple regions correspond to 68%, 95% CL regions drawn from the total event number in the excess 
region. The orange region is the one potentially compatible with the assumed total width. This 
bound is saturated when the partial widths from Oy, Oh are the only contribution to the total 
width. 


On the other hand, the hypothesis with only Xwz non-zero is highly disfavoured, with a 
Bayes factor of 2 • ICO 5 . 

We then proceed by drawing the best-fit regions for aww > G zz from the posterior 
L(a\yw,zz) 7T (^ww,zz)- The priors for the cross sections are taken flat. If one does not 
taken into the uncertainty on the background, the regions obtained are shown in the left 
pannel of Fig. 5. Finally, these regions of A ww,zz are readily translated into regions for 
total cross sections aww, azz, that are shown in Fig. 5. 

6 Interpreting the diboson excess with the neutral resonances EFT 

In the new physics scenarios considered in Sec. 3, the neutral resonance couplings to field 
strengths are universal. We will therefore make a simplifying assumptions and use a single 
parameter 

fv = fw = f b = fc i (6.1) 

both in the spin-0 and spin-2 cases. We then focus on the tree-level production induced 
by the Oh, Oy operators via gluon fusion (GGF) and weak boson fusion (VBF). We 
find that VBF is subleading to GGF for most of the parameter space. For the spin- 
0 case, these two operators also completely fix the width Tsm into SM particles (up to 
suppressed contributions from Or), while for the spin-2 case, one can have contributions 
from the operators O^, if present. i() In addition one can allow for an invisible width, 

16 In our analysis we do not take into account production of the spin-2 resonance via quark-fusion (which 
can be induced by the operators for ip = qL, Ur , dn), nor NLO-QCD effects. Both effects have recently 
been considered in Ref. [53] and were shown to lead to 0(1) modifications of the production rates. 
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i.e. Ttat = Tsm + Tiny into non-SM particles. The total width estimation from the shape 
analysis of Sec. 4.2 then provides us with another constraint in the fy — fu plane, i.e. 

Ty + T h < T t ot, (6.2) 

where the Th,v are the partial widths induced by the Oh, Oy operators. 

We write FeynRules [54] models for the EFT of neutral resonances described in Sec. 2. 
The signal expected from the spin 0 and 2 resonances pp —> (p —» WW,ZZ, and pp —» 
(f’/iu ~t WW, ZZ at the 8 TeV LHC are then computed using MadGraph 5 [55]. The main 
cuts are pt > 540 GeV and \r]\ <2 for each of the outgoing vector bosons, which we 
implement using MadAnalysis 5 [56]. 

The limits on the spin-0 CP even resonance are shown in Fig. 6. We choose T to t within 
the 95% confidence interval provided in Tab. 3, i.e. we fix it to T to t = 150 GeV (left panel 
of Fig. 6) and r to t = 20 GeV (right panel of Fig. 6). The orange shaded regions are those 
allowed by the condition (6.2). We also display how various scenarios fit into the shown 
parameter space. The Higgs portal scenario, which is indistinguishable from the radion 
in the RS brane model, is shown as the blue point. The line emerging from that point 
would correspond to a hypothetical model where the scalar boson can decay into invisible 
states. Similarly, the radion of the RS bulk model is depicted as the red point. Neither 
scenario can fit the observed excess, mainly because the operator Oq is not available for 
GGF production. Generating sufficient contribution from the VBF process would require 
too small values of fn, in conflict with the measured width. Finally, the green line shows a 
hypothetical scalar with universal fn = fv = /• This scenario could explain the required 
width and production rate for / = 7 TeV (19 TeV) for r tot = 150 GeV (20 GeV). 

We present the analogous limits on spin-2 resonances in Fig. 7. The RS brane model is 
shown as the green point corresponding to the value of k obtained from the chosen width 
of the resonance, again 150 GeV (left) and 20 GeV (right). The implicit values for the 
coupling are k = 0.23 (for r tot = 150 GeV) and k = 0.09 (for r to t = 20 GeV). Note that 
Th + Ty ~ 0.52 T t ot, the rest being contributed by the fermions. The RS bulk model with 
universal brane kinetic terms for the gauge fields has two parameters (rq and k) , and hence 
for fixed width it corresponds to one-dimensional curves in the parameter space, shown as 
the red curves. We vary 0 < rq < 90 (see footnote 10), corresponding to implicit values 
of k = 1.09 — 0.42 (for r tot = 150 GeV) and k = 0.40 — 0.16 (for r tot = 20 GeV). One 
can see that the brane model can in fact explain the observed excess, while the canonical 
RS bulk model (with rq = 0, the blue point) cannot, because the values of k needed to 
fit the correct width and production cross section are in conflict with each other. On the 
other hand, allowing for IR BKT’s, one can fit both width and total production rate of the 
excess. The required size of the BKT’s is only rq ~ 1 — 4, depending on r to t. 

7 Conclusion 

New particles, singlets under the SM interactions and with masses near the TeV scale can 
arise in many well-motivated extensions of the SM, including extra dimensions, strongly- 
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Figure 7. Constraints from the ATLAS diboson excess [7] on a generic neutral resonance with 


spin 2. A total width of 150 GeV and 20 GeV is assumed on the left and right pannel respectively. 


Purple regions correspond to 68%, 95% CL regions drawn from the total event number in the excess 
region. The orange region is the zone potentially compatible with the total width assumed. This 


bound is saturated when the partial widths from Oy, Oh are the only contribution to the total 
width. 


coupled scenarios as well as the Higgs portal. Such particles can be linearly coupled to SM 
operators, and can thus appear as resonances in s-channel processes. 

In this paper we first lay down the complete effective Lagrangians for neutral resonances 


of spin 0 and 2. It turns out that this EFT consists of only few operators, which can further 


be restricted by theoretically well-motivated assumptions, such as approximate flavor and 
CP conservation. Given the concise description of a large class of models in terms of few 
parameters, our EFT can serve as a model-independent framework to study the implications 
for any resonance searches at the LHC. We compute the generic widths of the resonances 
and present explicitly the matching to the new physics scenarios quoted above. 

We then investigate the possibility that a new heavy resonance be at the origin of the 
excess in diboson production recently reported by the ATLAS collaboration. We compute 
the local significances in both frequentist and Bayesian frameworks, showing a moderate 
evidence for the existence of a signal. We perform a shape analysis of the excess under full 
consideration of the systematic uncertainties to extract the width T to t of the hypothetical 
resonance, finding it to be in the range 26 GeV < Ttot < 144 GeV at 95% C.L. Turning 
to the study of total event numbers, we first evaluate the conditional probabilities for 
tagging a true W, Z and QCD jet as either W or Z from the ATLAS simulations. From 
these one deduces the tagging probabilities for the WW, WZ, ZZ selections reported by 
ATLAS. We further observe that these three overlapping selections follow a joint trivariate 
Poisson distribution, which opens the possibility of a thorough likelihood analysis of the 
event rates. The tagging probabilities are checked against the full observed sample. A 
conservative treatment of the dijet background is adopted, that includes the correlations 
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among the three selections WW, WZ, ZZ. The uncertainty on this background estimation 
is then taken into account as a systematic error in the total likelihood. Finally, using an 
actual hypothesis testing, we show that the data do not require WZ production and can 
thus in principle be explained by neutral resonances. 

Finally, we test the effective Lagrangians of neutral resonance using both the infor¬ 
mation from the width and the cross section of the analysis of the ATLAS data. It turns 
out that these pieces of information imply stringent contraints on the EFT parameter 
space once put together, even after including the background uncertainty. These exclusion 
bounds further imply that various popular scenarios appear to be totally incompatible with 
the ATLAS diboson excess. We find that neither scalars coupled via the Higgs-portal nor 
the RS radion can explain the ATLAS anomaly. The RS graviton with all matter on the IR 
brane can in principle fit the observed excess, while the RS model with matter propagating 
in the bulk requires the presence of IR brane kinentic terms for the gauge fields. 

As an outlook, we emphasize that it would be interesting to constrain the EFT for 
neutral resonances using other LHC searches. As the effective Lagrangians are rather 
predictive, powerful conclusions can be expected by combining the information from various 
channels. 
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A Polarization tensors for spin-2 fields. 

The polarization tensors for massive spin-2 fields satisfy 


= 0, , e^ = 0, 




= e 




(A.l) 


and their orthogonality and completeness relations read [36] 


*s s pcr tt per *s s' pLV rss' 


(A.2) 


where the projector n was given in Eq. (2.13). They can be written in terms of tensor 
products of spin-1 polarization vectors 
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(A.3) 

We note that the fields h™ u in Ref. [57] are not canonically normalized, and correspondingly 
their polarization tensors (propagators) are \/2 (2) times ours. 
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B Partial widths 


In Sec. 2 we gave the total widths of the various resonances in function of the effective 
couplings. The partial widths, if required, can easily be obtained from these formulae. For 
the field-strength couplings, one can use the decomposition 


3 1 _ 2 
72 “ + 72 " — 72 “ + 
Jw JB Jw 


_|_ ^ w 


fw Jb 


+ 


s 2 c 2 

a W 

fw Ib 


+ ?c 2 s 2 


l 

fw 



(B.l) 


where the four terms correspond to WW, ZZ, 77 and Z 7 decays respectively For the 
coupling fn one has simply f^ 2 = \ffj 2 + \ff[ 2 + \ff { 2 , corresponding to (longitudinal) 
WW, ZZ, and hh decays. Finally, for the partial widths of the fermions one can use 
decompositions such as 


Nq L 1 NtR , Nb R 
fqL ftR fb R 



+ 3 



(B.2) 


corresponding to the partial widths of the top and bottom quarks, and similarly for the 
other SM fermions. 
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